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This Is a summary of the work accomplished during the first half 
of 1985 (Including the month of July). 

1. The principal goal of the project Is to Investigate the 
feasibility of rain volume estimation over fixed and floating areas 
using rapid scan satellite data following a technique recently devel- 
oped with radar data, called the Area Time integral (ATI) technique. 

To accomplish this task, ca:e studies were continuously selected 
(during the first half of 1985 as well as during 1984) on the basis of 
existing radar and satellite data sets which match in space and time. 

The radar and rapid scan GOES satellite data were collected during 
the Cooperative Convective Precipitation Experiment (CCOPE) and North 
Dakota Cloud Modifi cation Project (NDCMP). Six multi cell clusters and 
cells were analyzed to the present time. The small mesoscale areas of 
Austin and Houze (1972) or the mesoscale cloud systems described by 
Super and Heimbach (1980) as cells bigger than a cumulus cloud but 
smaller than a typical mesoscale convective complex were considered 
"clusters". Three clusters were selected on each day, 12 June and 
2 July, as Illustrated In Table 1. The 12 June clusters occurred during 
the daytime, while the 2 July clusters during the nighttime. A total 
of 86 time steps of radar and 79 time steps of satellite Images were 
analyzed. There were approximately 12-min time Intervals between radar 
scans on the average. 

The radar data were on hand at the Institute of Atmospheric Sciences, 
South Dakota School of Mines and Technology. Magnetic tapes of digital 
satellite data were acquired by the Department of Atmospheric Science, 
Colorado State University, from the Bureau of Reclamation, the Univer- 
sity of Wisconsin through the National Environmental Satellite Data 
Information Service (NESDIS), and the GOES data archive at CSU. The 
satellite data were processed at the IRIS/DRSES facility located at 
CSU. The radar data were processed at SDSM&T. Details concerning data 
processing were described In the first and second semiannual reports. 

The trends of the er.ho areas and of the temperature thresholds 
picked to match the radar return area are displayed In Fig. 1A, B, C, 

D, E, and F for the six cells, respectively. A visual Inspection of 
these graphs emphasizes a general similarity of evolution and a tendency 
for satellite products to lag that of radar (Fig. 1A, IB, and IE). A 
two-cycle oscillation emphasizing the multicell character of the clusters 
Is demonstrated, except for cell E. The cirrus spissatus debris Is 
visible In two of the cells, IB and X2. A hint of the first growing 
period Is present In the satellite data trend (In Fig. 1A), primarily 
because the clusters Initial stage was overlaid by debris from previous 
activity. 

Figures 2A, B, C, D, E and F display the maximum radar reflectivity 
(Z), the temperature of Infrared maximum count (T MX ), and the maximum 
echo height (H) for reflectivity values >20 dBz on an Interpolated 







OiUGMNAL in'*' 

TF POOR QUALITY 



«T riK u 0 II *0 BIT fix (I • II MIS 


Fig. 1: Trends of radar echo area (A) and Infrared temperature matching radar echo 

area (Tmx) for the six analyzed clusters. 
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step-by-step basis for all six cells. Again, as In Fig. 1, the 
similarity of the cluster evolution, the 30 minute time lag between 
radar and satellite peaks, and the cirrus debris problem are evident. 

It Is Important to note that during the presence of cirrus debris, the 
minimum cloud top temperature remains lower than initially. 

A summary of radar and satellite products for the six analyzed 
clusters Is presented in Table 2. Columns (2) to (5) represent radar 
products, column (7) to (10) Independent satellite products, while 
column (11) to (13) satellite products matching radar suantities. 

Columns (11) to (13) were computed for the storm lifetime using histo- 
grams. Two of such histograms are presented in Attachment A, Fig. 4. 

A visual Inspection of both radar and satellite products (Table 2) 
emphasizes variable geometries (physical appearance) of the six consid- 
ered clusters. The 2 July clusters are more extended in the vertical 
than the 12 June as viewed by satellite. The radar products are roughly 
the same for both days. The MEH (column 3) average is only slightly 
higher (.196 km) for 2 July. The cluster with the largest MEH (X 2 , 

MEH = 14.4 km) did not generate the greatest amount of rain. 

The main purpose of this work is to reach to compute convective 
rain volumes over areas by application of the ATI technique based only 
on satellite data. Pie key element is to determine ATI from independent 
satellite data without consideration of radar reflectivities. This is 
possible if trends of satellite products generated independently (Table 2, 
columns 7-10) are similar to those of satellite based upon radar obser- 
vations as done here. For instance regression curves correlating satellite 
digital count thresholds to satellite IR maximum count values, or satellite 
IR maximum count values averaged over the duration of the cell have to be 
Identified. Such a linear regression relating the satellite IR maximum 
count value averaged over the duration of the cluster to the satellite IR 
temperature which matches the radar ATI is displayed in Fig. 5 of Attach- 
ment A. Certainly this finding is based only on a very limited sample 
size. Investigations are to further test and refine these results. 

Descriptive models representing the structure of the radar 
reflectivities and the minimum temperature of infrared digital counts 
as a function of surface area and time have a similar appearance above 
a given satellite threshold (Fig. 6A, B; Attachment A). This corres- 
pondence is emphasized by the -50°C and the 35 dBz curves. These 
models were computed by averaalng the six analyzed cluster products. 

For warmer thresholds, the muitlcell structure of the clusters gener- 
ates multimode digital count trends. The multicell structure of the 
cluster appears to be more evident on the satellite products. 

Adjustments of the radar and satellite data on a step-by-step 
basis were performed on 2-3 May 1985 during a joint SDSM&T and CSU 
meeting at Rapid City. Participants were T. H. Vonder Haar and 
P. Laybe from CSU and A. A. Doneaud, J. R. Miller, Jr., and L. R. 

Johnson from SDSM&T. 


Characteristic Products of Bowman Radar and GOES Satellite Data for SI* Clusters (or Cells) of Summer 1981 
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Some other 4-6 clusters are scheduled to be processed and analyzed 


next. 


2. A contribution entitled, "The Area-Time- Integral Technique 
to Estimate Convective Rain Volumes Over Areas Applied to Satellite 
Data — A Preliminary Investigation" coauthored by Andre A. Doneaud, 
James R. Miller, Jr., and L. Ronald Johnson of the Institute of Atmo- 
spheric Sciences, South Dakota School of Mines and Technology, Rapid 
City, South Dakota, and Thomas H. Vonder Haar and Patrick Laybe, of the 
Department of Atmospheric Sciences, Colorado State University, Fort 
Collins, Colorado, has been accepted for presentation at the Sixth 
Conference on Hydrometeorology of the American Meteorological Society 
29 October - 1 November 1985, Indianapolis, Indiana. A copy of the 
Preprint Volume Paper Is attached as Attachment A. 

3. A 15 minute presentation was given by A. A. Doneaud at the 
Panel Review Meeting of the Global Scale Atmospheric Research Program, 
8-11 July 1985, Columbia Inn, Columbia, Maryland. A trip report Is 
attached as Attachment B. An abstract of the presentation was 1r eluded 
In an "Abstract Digest" volume published by Science and Technolojy 
Corporation. The abstract Is attached as Attachment C. 

4. A 50 minute seminar was given by A. A. Doneaud at the NASA/ 
GSFC on the afternoon of 11 July. The subject was related to the 
progress made up to date on Grant No. NAG 5-386. 

5. A journal manuscript to be submitted for possible publication 
to the Monthly Weather Review or to the Journal of Cl imate and Applied 
Meteorology is under work and will be ready in 1-2 months. 
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THE AREA-TIME-INTEGRAL TECHNIQUE TO ESTIMATE CONVECTIVE RAIN VOLUMES OVER AREAS 
APPLIED TO SATELLITE DATA - A PRELIMINARY INVESTIGATION 


Andre A. Doneaud, James R. Miller, Jr., and L. Ronald Johnson 

Institute of Atmospheric Sciences 
South Dakota School of Mines and Technology 
Rapid City, South Dakota 57701-3995 

and 

Thomas H, Vonder Haar and Patrick Laybe 

Department of Atmospheric Sciences 
Colorado State University 
Fort Collins, Colorado 80523 


1. Introduction 

The main emphasis of this study is to 
investigate the f easibility of rain volume estima- 
tion using satellite data following a techniqut 
recently developed with radar data, called the 
Aea -Time-Integral (ATI) technique. The basis of 
the method is the existence of a strong correla- 
tion between the radar echo area coverage 
integrated over the lifetime of the storm and the 
radar estimated rain volume. It provides a means 
of estimating total rain volumes over fixed and 
floating target areas of the order of 1,000 to 
100,000 km 2 for clusters lasting >40 min (Doneaud 
et al . . 1981, 1984a). The technique does not 
require the consideration of detailed interior 
structure measured by radar or other means to 
generate rain volumes; only the area covered by 
rain events. Simplicity of the approach, if 
applied successfully to satellite data, holds the 
promise of reducing sources of error by minimizing 
required data input in comparison to other 
methods. 

Some of the satellite methods used to 
estimate rain volumes over areas considered cloud 
histories in the calculation process (i.e., 
Griffith et al . , 1976). Negri et al . (1984) found 
the Griffith-Woodley Technique TSwT) quite compli- 
cated as it necessitates different conversion 
factors for each stage of cloud development for a 
given climate. The cited authors suggested a 
simplified version of the technique. Griffith et 
al . (1981) showed that special precautions should 
be used to adjust the magnitude of the estimates 
when using the GWT technique in a much drier cli- 
mate. Scofield and Oliver (1977) have developed 
rainfall estimation methods for use in forecast 
offices primarily for flash flood prediction 
(large and long-lasting storms). Negri and Idler 
(1981) and Adler and Mach (1984) found cloud 
height-rainfall rate relations by considering 
only the period of maximum development of the 
convective system 


The ATI technique considers cloud h'stories 
and requires only one conversion factor to convert 
the ATI to total rain volume. The proposed method 
limits input to only satellite variables in esti- 
mating rain volumes. Relations between digital 
count thresholds which yield areas that match 
areas >25 d8z used to determine the ATI's havp to 
be found to extend the technique. Otner satellite 
quantities such as cloud top minimum temperature 
or thermodynamic profiles may also be useful to 
better understand selection processes that occur. 

2. Data and data processing 

Radar and rapid scan satellite data were 
collected during the summer season of 1981. The 
radar data were collected in Bowman (western North 
Dakota) as a part of the North Dakota Cloud 
Modification Project (NDCMP) and the rapid scan 
satellite data were collected as a part of the 
Cooperative Convective Precipitation Experiment 
(CC0PE). The radar data were processed at the 
South Dakota School of Mines and Technology and 
the satellite data at the Colorado State 
University. 

The Bowman radar was an Enterprise 
Electronics Corporation WR-100, 5.4-cm system 
equipped with a digital video integrator and 
processor and a magnetic tape recorder. The data 
tapes were processed following procedures similar 
to tho'-e described by Scliroeder and Klazura 
(1978). The rain volume for every cluster was 
computed using an optimized Z-R relationship, 

Z * 155 R* • 88 (Smith et^ al., 1975) based on data 
from this region. Oetai Is concerning the pro- 
cessing of the radar data can be found in Doneaud 
et al. (19d4c). Digital printouts of the dBz 
values at low tilt angle were prepared and con- 
verted from spherical to a rectangular coordinate 
system. 

Digital Visible and Infrared Spin Scan 
Radiometer (VISSR) data from both the east and 


west GOES were used. GOES-West was In the rapid 
scan (high frequency imaging) mode during the days 
analyzed. Rapid scan imagery from GOES-West was 
used to match as closv. as possible the satellite 
image times to the radar scan times. When a close 
match was not possible with GOES-dest imagery, 
GOES-East imagery was employed. The satellite 
data were processed at the IR1S/0RSES (Interactive 
Research Imaging System/01 rect Readout Satellite 
Earth Station) facility located at CSU. The 
satellite imagery was remapped into a Lambert con- 
formal projection so that the radar PPI data and 
the satellite imagery are in the same perspective. 
The resultant image resolution is -2.47 km 2 per 
image pixel. The visible satellite data was nor- 
malized for sola r zenith angle. The infrared 
satelluw data h.is a calibration stability of 
15° at 295°K. Tie remapped infrared imagery was 
modi fed by a routine which smoothes the sharp 
contrast between pixels by averaging. 

A radar sector of interest was defined tv 
delineate specific radar ecno clusters for each 
radar time throughout the radar echo cluster life- 
time. The radar sector of Interest was used to 
locate the convective cluster responsible for the 
rainfall in the remapped satellite data. The 
radar echo clusters within the defined sectors of 
interest were identified by drawing a “box" using 
axes and radials around the cluster. The cluster 
echo areas of >25 d&z reflectivity threshold and 
the corresponding radar AT! and ra i r, volur.e are 
calculated. 

Next, a satellite sector of interest ias 
defined by applying small adjustments to the radar 
sector using a manual processing technique. This 
was done to avoid cloud features suspected of not 
being detected by the radar. It is one of the 
most delicate tasks of this investigation because 
of several considerations: a) satellite and radar 

systems respond to different chararteristics of 
clouds at different atmospheric levels; b) time 
differences betveen satellite and radar data sets 
can vary by as '■och as 10 min; c) spatial posi- 
tioning of the satellite observations relative to 
the radar location necessitate error inclusion due 


to limited accuracy of geometric corrections; 
d) vertical wind shear advects the cloud tops 
downwind from the location of the radar echo area; 
and e) in the definition of the matching satellite 
sector of interest, tne inclusion or exclusion of 
cirrus spissatus debris is a very important deci- 
sion. If the cirrus debris appeared to be comple- 
tely detached horizontally from the radar echo 
cluster, the cirrus cloud cover was excluded from 
the satellite sector of interest. 

Six clusters lasting 2-4 hours are analyzed 
for which both digital satellite and radar data 
existed. These storm events occurred on 12 June 
and 2 July. Table 1 describes some characteristic 
quantities of these clusters. The 12 June 
dusters occurred during the daytime, while the 
2 July dusters occurred during the nighttime. A 
total of 36 time steps of radar and 79 time steps 
of satellite images were analyzed. There were 
approximately 12-min time Intervals between scans 
on the average. The 12 June clusters were ana- 
lyzed using both Infrared and visible data. For 
the 2 July clusters, only infrared data were 
available. The small mesoscale areas of Austin 
and Houze (1972), or the mesoscale cloud systems 
described by Super and Heimbach (1980) as cells 
bigger than a cumulus cloud but smaller than a 
typical mesoscale convective complex were 
considered “clusters" in this work. 


3. Analysis of the data 

3.1 Time Step Analysis 

Using the technique described in Section 2. 
satellite brightness and temperature pixel counts 
above defined thresholds within the satellite 
sector of interest were counted and converted to 
area. Histograms of area versus brightness counts 
were produced for each satellite time step for 
each data type for the lifetime of the storm. 

These histograms and the computed radar area 
(multiplied by the radar time interval) were used 
to determine thresholds for satellite data 
(matching the radar area multipled by time) by 


TABLE 1 

Characteristics of the Analyzed Clusters 


No. 

Cluster 
Ident . 

Date 

1981 

Time Ppriod (Radar) 

Time 

Radar 

Steps 

Satellite 

Cluster Type 

Satell ite 
Data Type 




(GMT) 





1 

1A 

12 June 

1657 - 1955 

14 

13 

Multi cel 1 

Infrared, Visible 

2 

IB 

12 June 

1746 - 1955 

1G 

12 

Multi cel 1 

Infrared, Visible 

3 

E 

12 June 

1657 - 1930 

11 

12 

One Cell 

Infrared, Visible 

4 

XI 

2 July 

0103 - 0453 

19 

16 

Multicell 

Infrared 

5 

X2 

2 July 

0201 - 0434 

13 

10 

Multicell 

Infrared 

6 

A 

2 July 

0103 - 0453 

19 

16 

One Cel 1 

Infrared 
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forcing the areas multiplied by time for satellite 
to be equal to radar. Both satellite and radar 
time series were reduced to a constant 15-min time 
increment by linear interpo 1 at ion such that func- 
tional time deoendence is comparable. The echo 
areas >25 dBz (multiplied by time), as well as the 
corresponding threshold count values, exhibit a 
large spectrum of variations during the storm 
lifetime. No relationship has been found between 
these variations and the stages of the storm 
lifetime. 

The trends of the echo areas and of the 
temperature thresholds picked to match the radar 
return area are displayed in Fig. 1A, B, C, 0, E, 
and F for the six cells, respectively. A visual 
Inspection of these graphs emphasizes a general 
similarity of evolution and a tendency for 
satellite products to lag that of radar (Fig. 1A, 
IB, and IE). A two-cycle oscillation emphasizing 
the multiceli character of the clusters is 


demonstrated, except for Cell E. "he c'rrus 
spissatus debris is visible in two of the cells, 

IB and X2. fhe cirrus debris maxe the definition 
of the end of a cluster (as a convective entity) 
from the satellite Images difficult. A h'nt of 
the first growing period 's present in the 
satellite data trend (in Fig. 1A), primarily 
because the clusters initial stage was overlaid 
by debris from previous activity. Again remnants 
of previous activity preclude initial observa- 
tions in Fig. 10 (Cluster XI) and Fig. *F 
(Cell A). 

Figure 2A. B, C, 0, i, and F dually the 
maximum radar reflectivity ( Z ' , the temperature of 
infrared maximum count (Tvix), and tne maximum echo 
height (H) for reflectivity values >20 dBz un an 
interpolated step-by-step basis for all s*< cells. 
Again, as in Fig. 1, the similarity of cluster 
evolution and the time lag between radar and 
satellite is evident. 
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s- < analyzed clusters. 
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Trends of radar maximum echo height (H), maximum reflectivity (Z). and temperature of infrared 
count (Tmx) for the s ’ x ® ni 'y Je ' 1 Clusters. 


In Fig. 1A, a 30-min lag between the 
maximum echo height and maximum reflectivity as 
radar products and the minimum temperature by 
satellite is evident. That again demonstrates 
the multicell character of this cluster ( 1 A ) . An 
identical 30-min lag is also visible in Figs. 28 
and 2C. In tie same two figures, the cirrus 
debris is evident. It is Important to note that 
during the presence of cirrus debris, tho minimum 
cloud top temperature remains lower than ini- 
tially. Remnants of previous activity precluded 
definitive in 1 1 a 1 observations for Cells XI and 
X2 (Figs. 20 and 2E). For Cell X2, only the 
infrared maximum count tr«*nd converted to 
temperature is available. 

The step-by-step ladar reflectivites for 
>25 dU* on a 10-dBz step interval for the 2 July 
1881 clusters, as well as the corresponding 
satellite IR features, are present'd >n Fig. 3. 
The radar sectors used to delineate the echo 


clusters are displayed on each radar PPI image, 
while the satellite subjective sectors are 
displayed on the satellite -mages. The IH .on- 
tours represent -4°C temperature steps with the 
first contour representing a fixed tenperatu. e. 
The multicellular structure of Cluster XI is evi- 
dent. Cluster X2 starts as a single cell and 
develops into a multAel lular structure as it 
matures. In the satellite Images as well as in 
the radar PPI, Cell <1 grows until 0128 GMT (two 
consecutive time steps); it slowly decays until 
0241 GMT, and nure rapidly until the end at 
0503 GMT. Cell A is a part of Cluster XI, which 
grew until 02<1 GMT, then slowly decayed two con- 
secutive time steps. The evolution of Cell <2 is 
best followed by the radar PPI. The data suggest 
that the development of cell X2 is *n the outflow 
region of a previous cluster, well »xtended hori- 
zontally and vertically. The maximum vertical 
extension of all clusters occurs during the 
first half of their lifetime. 
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Cluster LI fetime Analysis 

A summary of radar and satellite products 
for the six analyzed clusters is presented in 
Table 2. 

Columns (2) to (5) represent radar 
products, while columns (7) to (13) are satellite 
products. The digital count thresholds matching 
radar ATI's, and the respective minimum tem- 
peratures, as well as the satellite ATI (columns 
11-13) were computed for the storm lifetime using 
histograms. Two of the histograms are presented 
for Cluster 1A (VIS and Id); Fig. 4A, B. The 
dashed-dotted line represents the satellite deter- 
mined ATI (on a step-by-step basis) and its 
respective digital count. The columns on Table 2 
displaying satellite products could be divided in 
columns representing independent satellite pro- 
ducts [columns (7), (8), (9), (10)], and columns 
representing satellite products obtained as a 
combined radar versus satellite analysis [columns 
(11), (12), (13)]. 

A visual inspection of both radar and 
satellite products (in Table 2) emphasizes 
variable geometries (physical appearance) of the 
six considered clusters. The 2 July clusters are 
more extended in the vertical than on 12 June as 
viewed by satellite. The satellite products 
(columns 7-12) are always greater for 2 July than 
12 June. The radar products are roughly the same 
for both days. The MEH (column 3) average is only 
slightly h'gher (0.96 km) for 2 July. The cluster 
with the largest MEH (X2, MEH = 14.4 km) did not 
generate the greatest amount of rain. 

The six clusters exhibit a large spectrum 
of ATI values ranging from 704 km^ hr for Cell A 
to 4,625 km^ hr for cluster 1A. The rain volume 
varies from 4,082 km 2 m (IB) to 16,119 km^ mm 
(1A). A large diversity in cluster size is not 
demonstrated by this data set. 


The main purpose of this investigation is 
to compute convective rain volume over a moving 
targeted area by application of the ATI technique 
based only on satellite data. As such, the key 
element is to determine ATI from independent 
satellite data without consideration of radar 
data. This is possible if trends of satellite 
products (Table 2, columns 7, 8, 9, and 10) 
generated independently are similar to those of 
satellite products based upon radar observatons 
as done here. An independent satellite quantity 
which defines the most appropriate digital count 
that matches the radar ATI for a given cluster 
would be a solution. For instance, regression 
curves correlating satellite digital count 
thresholds to satellite IB maximum count values, 
or satellite IB maximum count values averaged 
over the duration of the cell have to be iden- 
tified. Such a linear regression relating the 
satellite >R maximum count value averaged over 
the duration of the cluster to the satellite IB 
temper, ture which matches the radar ATI is 
displayed in Fig. 5. The correlation coefficient 
is .93. Certainly tnis finding is based only on 
a very limited sample size. Investigations are 
anticipated to further test and refine these 
findings. 

Encouraging ; s the fact that descriptive 
models representing the structure of the radar 
reflectivities and the minimum temperature of 
infrared digital counts as a function of surface 
area and time have a similar appearance above a 
given satellite threshold (Fig. 6A,3). This 
correspondence is emphasized by the -5Q°C and 
the 35-dBz curves. These models were computed 
by averaging the six analyzed cluster products. 
For warmer thresholds, the multicell structure 
nf the clusters generates multimode digital count 
trends. The multicell structure of the cluster 
appears to be more evident in the satellite 
products. 


TABLE ! 

Chirictarl ttlc Product} of Bonn «M*r «»« SOES SltolHtt OjtA for Sll Clultart (or Colli) of Suiwor 1981 


■ADA» QATA SATELLITE 3AfA 


c#n io. 

And for at ion 
(hn) 

*ai IB: 

Echo 

km 

ATI 

hr 

Rain Vol 
km J 

A 

*4«1«U« 

Distil Count 

T*-so Of 
IR *4t Count 
T« (-•C) 

Avinw* 

T^> of !R 
*41 Count 
T AH! (-•Cl 

4i*r|fjod* 
VM *4i Count 

Digital 

Count Thrath 
*atcning Radar 
ATI 

Te-U *4tcn*nq 
Radar ATI 

T ma -’Cl 

*Tl 

(Sita.lital 

— 

(1) 

(2) 

(3) 

(4) 

(S) 

(8) 

(?) 

(8) 

(9) 

(101 

HD 

(12) 

(13) 

1* 

m 

61 

12.2 

4,675 

16.119 

tR 

VIS 

195 

66 

50.2 

41.5 

62.8 

179.70 

57.73 

33.0 

4.534 

4.534 

IB 

r. to 

V 

12.1 

1,11? 

4,082 

IR 

-IS 

194 

65 

49.2 

•5.1 

63.0 

189.09 

59.36 

45.0 

1.169 

1.169 

c 

2.T5 

61 

12.2 

1,302 

6,831 

IR 

VIS 

195 

64 

50.2 

43.6 

61.0 

185.68 

57.85 

41.0 

1,221 

1.221 

il 

03 

57 

12.5 

1.897 

11.001 

IR 

204 

59.2 

53.5 


201.21 

56.4 

1,376 

I? 

2775 

55 

14,4 

1.338 

1.027 

IR 

207 

61.2 

58.6 


202. 

58.0 

1,243 

A 

3733 

S3 

12.5 

704 

4 268 

IR 

2f.4 

59.2 

53.8 


201.39 

56.5 

>2 2 


"A.arjg, tit coll llfft’oo. 



ORIGINAL PSitz. ,c 
OF POOR QJALHY 




Kiq, 4 : Histograms of VIS (A) and IP ( B ) digital counts vs. satellite ATI's for cell 1A; the radar ATI 

and tne correspondent digital count are also displayed. 



fliLl: Scatter and linear regression of 

temperatures matching radar ATI and average 
(over the cell duration) temperatures of IP 
maximum count. Every stai represents a 
cluster. 


4. Concluding remarks 

This analysis of six convective clusters 
suggests that rhe extension of the ATI technique 
to the use of satellite data holds promise. The 
differences of the internal structure of the radar 
reflectivity features (taken close to the cloud 
base) and of the satellite features (taken near 
the top of the cloud) g ; ve rise to differences 
in estimating rain volumes by delineating area; 
however, by focusing uoon the area integrated 
over the lifetime of the storm, tne investigations 
suggest that some of the errors produced by the 
differences in the cloud geometries as viewed 
by radar or satellite a^e eliminated. 

Models describing the structure of radar 
reflectivities and minimum temperatures by 
infrared digital counts as a function of surface 
area and time are similar above given satellite 
and radar thresholds. 9y averaging the six 
analyzed cluster products, these thresholds were 
found to be -50°C and 35 dBz for satellite and 
radar data, respectively. 


radar reflectivity structure 


SATELLITE OIIITAL COUNT STRUCTURE 



Fig. 6: Structures of radar reflectivities (A) and minimum temperature of IH digital counts (8) as 

functions of normalized area an! time. 
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It is of crucial importance in this type 
of analysis to be able to determine accurately the 
start and the ena times of the convective events. 
The starting time can easily be determined when 
the first cluster's echo is visible in the radar 
data. Yet, the cirrus debris makes the definition 
of the end of the cluster (as a convective entity) 
from the satellite images difficult. This incon- 
venience can be overcome if the maximum develop- 
ment of the cluster is taken as the ending time of 
the analysis. The ATI technique also works if 
only the growing part of the cluster lifetime is 
considered (Doneaud et al . . 1984b). The maximum 
development of the cluster can be defined using 
radar data as the maximum echo height, the maximum 
reflectivity, or the echo area >25 dBz. Satellite 
data definition of maximum development could be 
based on tne maximum count value or the minimum 
temperature. Investigations are anticipated to 
further illuminate such a procedure. 

In some of the graphs presented (Figs. 1 
and 2), time lags as great as 30 min were 
observed. Maximum echo height, maximum reflec- 
tivity, and maximum echo area were found to occur 
before the IR temperature threshold matching echo 
areas >25 dBz and the minimum temperature of 
infrared maximum count. These time lags do not 
generate problems in considering the maximum 
development of the cluster as the ending time 
of the analysis. 

Samples of radar and satellite products 
averaged over the six analyzed clusters were 
compared. A linear regression relating the 
satellite IR maximum count value (as an indepen- 
dent satellite product) to the satellite IR 
temperature matching the radar ATI (dependent 
on radar - ) showed a high correlation (.93). This 
would indicate that definition of the ATI's value 
could be based on independent satellite data. 

The regression is based only on a very limited 
sample size, making concluding remarks hazardous. 
Investigations are anticipated to further test 
and develop these findings. 
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ATTACHMENT B 


31 July 1985 


MEMORANDUM FOR THE RECORD 
FROM: A. A. Dorfe^ 

SUBJECT: Trip Report to Attend the Panel Review Meeting of the NASA 

Global Scale Atmospheric Research Program, Washington, DC 


I traveled to Washington, DC, by private auto between 5-20 July 
1985 to attend the panel review meeting of the Global Scale Atmospheric 
Research Program. The panel review meeting was scheduled to be held 
at the Columbia Inn, Columbia, Maryland. The logistics of this review 
were handled by Science and Technology Corporation (STC) under the 
chairmanship of Dr. Karl Johannessen, STC. 

In addition to the participation at the panel review meeting I 
gave a seminar at the NASA Goddard Space Flight Center on the afternoon 
of 11 July and I visited NSF offices the afternoon of 9 July and the 
morning of 12 July. 

The panel review was conducted In an open forum similar to an AMS 
topical meeting. Since the Global Scale Atmospheric Research Program 
covers 60 tasks and each task was Initiated by an individually identi- 
fiable project, presentations were not combined and each presenter had 
15 minutes to show his objectives, results and future worx. 

The panel members were Dr. Richard Anthes (NCAR), Dr. David Burridge 
(ECMWF), Dr. John Hovermale (NEPRF), Dr. Karl Johannessen (STC, Panel 
Moderator), Pr. Donald Johnson (University of Wisconsin), Dr. Peter 
Webster (Penn State University), and Dr. Axel Winn-Nielsen (Danske 
Meteorologiske Institut). The Moderator's goal was to let the panel 
look at the progress made In each individual task and assess how that 
progress contributes to the overall balance of the program. In such a 
way redundant or nonproductive research could be removed and new tasks 
in research areas more beneficial to the program could be initiated. 

The technical program for the panel review was organized in three 
sessions: Session A-Satellite Data Retrieval and Objective Analyses, 
Chairperson K. R. Johannessen, ETC; Session B-General Circulation, 
Modeling and Global Diagnostic, Chairpersons E. Kalnay, NASA/GSFC and 
G. H. Flchtl, NASA/MSFC; Session C- Theoretical Studies with Special 
Models, Regional Diagnostic Studies, Teleconnections, Chairpersons 
R. K. Kakar, Jet Propulsion Laboratory and S. H. Melfi, NASA/GSFC. 

An opening session was chaired by Dr. K. R. Johannessen. Dr. J. S. 
Theon presented the NASA Headquarters viewpoint and Dr. Robert J. Curran 
(The Program Manager) presented the Program (of the Global Scale 
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Atmospheric Processes Research) Orientation. Dr. Johannessen Introduced 
the panel members. Following E. Kalnay (NASA/GSFC), G. Fichtl (NASA/MSFC), 
and R. K. Kakar (NASA/JPL) presented overviews of their institution 
program. 

Session A was dominated by papers dealing with retrievals of wind, 
temperature, moisture, and pressure fields. Most of the presentations 
were given by NASA/GSFC scientists. The GLAS (Goddard) retrieval model 
uses three major components: analysls/fcrecast/retrieval . Several 
significant developments, including the use of "significant height" 
orography, higher resolution and improved parameterization of super- 
saturation precipitation have been implemented. The GLAS Analysis 
System (Baker, 1983) has also undergone significant improvements, 
mainly the replacement of the two-dimensional by a three-dimensional 
scheme. The Goddard Temperature Retrieval System (Susskind et al . , 

1984) determines atmospheric and surface parameters which, wTien 
substituted in the radiative transfer equations, mar.ch the infrared 
and microwave observations within a specified amount. This physically 
based retrieval system determines, together with atmospheric temperature 
soundings several other weather and climate parameters. Significant 
recent improvements include the development of accurate sea surface 
temperature determination (Susskind and Reuter, 1985) and of humidity 
soundings. The goal of this Goddard Research Group Projects is to 
include in addition to the customary atmospheric fields available from 
the four-dimensional analysis (pressure, temperature, humidity, and 
winds) other fields like estimates of cloud cover, cloud top pressure 
and temperature, atmospheric heating rates, snow and ice cover, sur- 
face fluxes of heat, and others. Besides Goddard scientists, some 
University scientists presented their progress. G. L. Achtemeier et 
al. (Illinois) presented Variational Objective Analyses for Cyclone 
Studies; G. Wahba and D. R. Johnson (Wlsconsln-Madison) described some 
Mathematical Methods for Variational Objective Analysis; and M. J. 

Munteanu (Maryland) showed applications of Clustering Techniques and 
Total Ozone Data to Satellite Soundings. Our presentation. The ATI 
Technique (To Estimate Convective Rain Volumes Over Areas) Applied to 
Satellite Data, was included in Session A and was well received. 

Session B papers dealt mostly with accomplishments on general 
circulation modeling. Work on NASA/GSFC has been done mostly to 
improve and add new data to the fourth order G.C.M. available at NASA. 
Numerical diagnostic studies of major weather events were developed by 
R. Atlas. E. Kalnay et^al. explored the existence of atmospheric 
predictability in the montFTly to seasonal time scale by means of the 
use of comprehensive atmospheric models. The authors adapted the Lagged 
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Average Forecasting Techniques (LAFT) of Hoffman and Kalnay (1983, 1984) 
to 10 day operational forecasts. Y. Mlntz et aK_ designed a biophysical 
numerical model of the vegetated surface of~the earth which can be used 
with atmospheric general circulation models. H. M. Helfand developed 
the 6LAS fourth order forecast model (Kalnay et al., 1983). The devel- 
opment included improvements to the model's horizontal and vertical 
resolution, design of better high-latitude and Shapiro filters, intro- 
duction of a nonlinear normal modes initialization scheme, enhancement 
of the model's orographic boundary data sets and creation of more sophis- 
ticated schemes for the parameterization of sub-grid scale physical 
processes. In addition to NASA scientists, some university people 
presented their results. S. E. Cohn et al . (Courant Institute of Mathe- 
matical Sciences) developed an ImplicTt scheme for global numerical 
weather prediction. The method uses a large time step. Current primi- 
tive-equation models based on explicit and semi-implicit methods, must 
use a small time step to ensure numerical stability. D. A. Salstein 
and R. D. Rosen from Atmospheric and Environmental Research, Inc., 
Cambridge, Massachusetts, as well as Tsing-Chang Chen from Iowa State 
University determined physical processes generating the distribution 
of dlabatic heating throughout the atmosphere. Knowledge of the dis- 
tribution of dlabatic heating throughout the atmosphere is vital to 
understanding the way In which the energy cycle of the atmosphere is 
maintained through the generation of Available Potential Energy (APE). 
Tsing-Chang Chen used the thermodynamic equation in the isobaric coor- 
dinates to estimate the diabatic heating. His results generally agree 
with other previous investigations, even though some differences occurred. 

Session C was dominated by papers given by University scientists. 

B. Saltzmass (Yale University) developed models with nonlinear moist 
barotroplc and barocolinic processes Including viscosity and heating, 
and lower boundary orographic and thermal forcing. The authors related 
the dynamical activity In the atmosphere to satellite-derived infrared 
and visible radiation measurements in order to study atmospheric phe- 
nomena (l.e. hydrometeors). E. R. Reiter (Colorado State University) 
studied surface energy fluxes In complex terrain. J. A. Dutton et aU_ 
(Penn State University) studied global atmospheric dynamics in orcTer to 
develop inverse methods that use ctynamical criteria to convert satellite 
observations of phenomenological structure into quantitative observations. 
R. S. Lindzen (Massachusetts Institute of Technology) made a brief 
review of the studies carried out at MIT, studies of stationary waves, 
teleconnections, blocking, cumulus momentum transport, finite amplitude 
evolution of barocllnlc Instability, water vapor convergence due to 
temperature gradients and atmospheric free Rossby wave. J. P. McGuirk 
and A. H. Thompson (Texas A&M University) developed and Implemented a 
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statistical procedure for using TIROS N microwave date to Infer Infra- 
red channel data for overcast conditions. The authors found that the 
satellite-sensed temperature and moisture signals vary differently, 
but consistently, within moisture bursts. In an empirical orthogonal 
function analysis over the tropical eastern Pacific, they found that 
a) the vertically oriented eigenfunctions could be interpreted in terms 
of typical meteorological events while b) the horizontal distribution 
of the eigenfunction amplitudes relate these meteorological signals to 
moisture bursts. D. 6. Vincent (Purdue University) presented atmospheric 
circulation features In the southern hemisphere for the period 10-27 
January 1979 using Global Weather Experiment (FGGE) data. The study 
focused on one of the major features in that period, the South Pacific 
Convergence Zone (SPCZ). K. B. Katsaros (University of Washington) 
described the structure of midlatitude cyclones observed with the atmo- 
spheric water channels of the scanning multichannel microwave radiometer. 
S. Mudrlck (University of Missourl-Columbia) tested the quasi-geostrophic 
(QG) dynamics versus the primitive equation (PE) dynamics for modeling 
the effect of cyclone wave. By simulating a polar low possessing a 
wave length of 1 200 km the authors found a quick development of that 
low the first day (with a maximum of the eddy kinetic energy) and then 
a decrease to a relative minimum kinetic energy on the third day. 

The aforementioned papers were randomly selected without any 
technical or special considerations. A preprint voiume containing all 
of the abstracts of the papers is available for consultation. 

On the afternoon of 11 July I gave a seminar entitled "The ATI 
Technique and its Possible Application to Satellite Data" at the NASA/ 
GSFC. The seminar was organized by Andrew Negri from the NASA/GSFC 
Laboratory for Atmospheric Science. Around 40 participants attended 
the seminar which was well received by the attendees. 

I returned to Rapid City Saturday, 20 July, after taking three 
days of vacation. 
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ATTACHMENT C 


THE AREA-TIMt-INTEGRAL TECHNIQUE (TO ESTIMATE CONVECTIVE 
RAIN VOLUMES OVER AREAS) APPLIED TO SATELLITE DATA 

A. A. Doneaun, J. R. Miller, Jr., and L. R. Johnson 

Institute of Atmospheric Sciences 
South Dakota School of Mines and Technology 
Rapid City, South Dakota 57701-3995 

and 

T. Vonder Haar and P. Laybe 

Department of Atmospheric Sciences 
Colorado State University 
Ft. Collins, Colorado 35023 


EXTENDED ABSTRACT 

The main emphasis of the study is to investigate the feasibility of rain 
volume estimation using satellite data following a technique recently developed 
with radar data (collected in the northern High Plains), called the Area Time 
Integral (ATI) technique (Doneaud et al , , 1981; 1984a; 1984b). The technique 
has been verified and proven to be efficient in Florida (Lopez, 1983) and South 
Africa. The basis of the method is the existence of a strong correlation 
between the area coverage (>25 dBz) integrated over the lifetime of the storm 
(ATI) and the rain volume. The technique does not require the consideration of 
the structure of the radar ntensities las, for instance, the Z-R technique) to 
generate rain volume. This fact simplifies and may reduce some sources of error 
in applying the technique to satellite data. Similarly, it is suggested that 
once radiance thresholds matching ATI's are determined by using visible (VIS) 
and infrared (IR) information integrated over a storm's lifetime, the digital 
unit threshold should exhibit variation among storms only in relation to the 
geometry (i.e., vertical extension) of the convective cloud system; adjustments 
may be necessary for regional and seasonal differences. The ATI technique may 
also be used if only the growth part of the cluster is considered (Doneaud et 
al . , 1984b). That eliminates the difficulties related to the definition of the 
end of the cluster's lifetime (as a result of cirrus debris) using satellite 
data. The new technique, if proven valid for satellite data, may be cate- 
gorized as a simplified Woodley-Griffith technique adjusted for measurable 
characteristics of cloud systems (Negri and Adler, 1981). 

Four multicell clusters were analyzed in 1Q84. Cells bigger than a cumulus 
cloud, but smaller than a typical mesoscale convective complex, were considered 
"clusters." The radar data were processed at the South Dakota School of Mines 
and Technology. The digital rapid scan satellite data were processed at the 
I R I S/DRSES facility located at Colorado State University. Different radar and 
satellite quantities for every time step for the selected clusters were computed. 
They reflect the multicell characteristic of the clusters, and yet still exhibit 
similar trends. Some preliminary results were presented at the 1984 Conference 
on Radar Meteorology (Dcneaud et a 1 . . 1985). Six additional cells will be 
processed and analyzed in 1985. 


The key element of this investigation is to determine an ATI from satellite 
data without using radar reflectivities. A satellite quantity defining the most 
appropriate digital count matching the radar ATI for a given cluster needs to be 
found: for instance, correlation studies between satellite digital count 

thresholds (matching radar ATI's), and satellite IR maximum cloud values (°C) 
are planned. Thp clusters analyzed show that the digital counts matching the 
radar ATI's are related to cluster characteristics and environment conditions. 
Twelve to fourteen clusters are scheduled to be analyzed through this grant 
(ending in December 1986). 

A new proposal is mienuc.. be submitted for the period 1987-89 to verify 
the applicability of the technique in regions with different climatic conditions 
using independent sets of data (i.e., data from the Oklahoma "Pre-STORM" 
Experiment). The objective of the new proposal should be twofold: to extend 

the validity of the technique in another region with possible extension to a 
global atmospheric scale, and to determine its accuracy using an independent 
data set. Sensitivity and accuracy tests will be proposed.- The new cases 
will be finally added to the total number of cases describing the technique. 


REFERENCES 

Doneaud, A. A., P. L. Smith, A. S. Dennis and S. Sengupta, 1981: A simple 

method tor estimating convective rain volume over an area. Water 
Resources Res. , 17 , 6, 1676-1682. 

, S. Ionescu -Ni scov, D. L. Priegnitz and P. L. Smith, 1984a: The 

area-time integral as an indicator for convective rain volumes. J . 

Climate App 1 . Meteor. , 23, 555-561. 

, , and J. R. Miller, Jr., 1984b: Convective rain rates 

and their evolution during storms in a semi-arid climate. Mon. Wea. 

Rev. , 112 , 1602-1612. 

, J. R. Miller, Jr., L. R. Johnson, T. H. Vonder Haar and P. Laybe, 

1985: An attempt to extend the ATI technique to estimate convective 

rain volumes using satellite data. Preprints 22nd Conf. Radar Meteor. , 
10-13 Sep 1985, Zurich, Switzerland, --. 

Lopez, E. R., J. Thomas, D. 0. Blanchard and R. L. Holle, 1983: Estimation 

of rainfall over an extended region using only measurements of the area 
covered by radar echoes. Preprints 21st Conf. Radar Me teor. , Amer. 

Meteor. Soc., Edmonton, Alberta, Canada, 681-681T 

Negri, A. J., and R. F. Adler, 1981: Relation of satellite based thunderstorm 

intensity to radar-estimated rainfall. J . App 1 . Meteo r. , 20 , 288-300. 


